Systemic interindividual epigenetic variation in humans is associated with transposable elements and under strong genetic control.

BACKGROUND: Genetic variants can modulate phenotypic outcomes via epigenetic intermediates, for example at methylation quantitative trait loci (mQTL). We present the first large-scale assessment of mQTL at human genomic regions selected for interindividual variation in CpG methylation, which we call correlated regions of systemic interindividual variation (CoRSIVs). These can be assayed in blood DNA and do not reflect interindividual variation in cellular composition. RESULTS: We use target-capture bisulfite sequencing to assess DNA methylation at 4086 CoRSIVs in multiple tissues from each of 188 donors in the NIH Gene-Tissue Expression (GTEx) program. At CoRSIVs, DNA methylation in peripheral blood correlates with methylation and gene expression in internal organs. We also discover unprecedented mQTL at these regions. Genetic influences on CoRSIV methylation are extremely strong (median R2=0.76), cumulatively comprising over 70-fold more human mQTL than detected in the most powerful previous study. Moreover, mQTL beta coefficients at CoRSIVs are highly skewed (i.e., the major allele predicts higher methylation). Both surprising findings are independently validated in a cohort of 47 non-GTEx individuals. Genomic regions flanking CoRSIVs show long-range enrichments for LINE-1 and LTR transposable elements; the skewed beta coefficients may therefore reflect evolutionary selection of genetic variants that promote their methylation and silencing. Analyses of GWAS summary statistics show that mQTL polymorphisms at CoRSIVs are associated with metabolic and other classes of disease. CONCLUSIONS: A focus on systemic interindividual epigenetic variants, clearly enhanced in mQTL content, should likewise benefit studies attempting to link human epigenetic variation to the risk of disease.

DNA methylation at a nutritionally sensitive region of the PAX8 gene is associated with thyroid volume and function in Gambian children.

PAX8 is a key thyroid transcription factor implicated in thyroid gland differentiation and function, and PAX8 gene methylation is reported to be sensitive to the periconceptional environment. Using a novel recall-by-epigenotype study in Gambian children, we found that PAX8 hypomethylation at age 2 years is associated with a 21% increase in thyroid volume and an increase in free thyroxine (T4) at 5 to 8 years, the latter equivalent to 8.4% of the normal range. Free T4 was associated with a decrease in DXA-derived body fat and bone mineral density. Furthermore, offspring PAX8 methylation was associated with periconceptional maternal nutrition, and methylation variability was influenced by genotype, suggesting that sensitivity to environmental exposures may be under partial genetic control. Together, our results demonstrate a possible link between early environment, PAX8 gene methylation and thyroid gland development and function, with potential implications for early embryonic programming of thyroid-related health and disease.

Identification of cell type-specific methylation signals in bulk whole genome bisulfite sequencing data.

BACKGROUND: The traditional approach to studying the epigenetic mechanism CpG methylation in tissue samples is to identify regions of concordant differential methylation spanning multiple CpG sites (differentially methylated regions). Variation limited to single or small numbers of CpGs has been assumed to reflect stochastic processes. To test this, we developed software, Cluster-Based analysis of CpG methylation (CluBCpG), and explored variation in read-level CpG methylation patterns in whole genome bisulfite sequencing data. RESULTS: Analysis of both human and mouse whole genome bisulfite sequencing datasets reveals read-level signatures associated with cell type and cell type-specific biological processes. These signatures, which are mostly orthogonal to classical differentially methylated regions, are enriched at cell type-specific enhancers and allow estimation of proportional cell composition in synthetic mixtures and improved prediction of gene expression. In tandem, we developed a machine learning algorithm, Precise Read-Level Imputation of Methylation (PReLIM), to increase coverage of existing whole genome bisulfite sequencing datasets by imputing CpG methylation states on individual sequencing reads. PReLIM both improves CluBCpG coverage and performance and enables identification of novel differentially methylated regions, which we independently validate. CONCLUSIONS: Our data indicate that, rather than stochastic variation, read-level CpG methylation patterns in tissue whole genome bisulfite sequencing libraries reflect cell type. Accordingly, these new computational tools should lead to an improved understanding of epigenetic regulation by DNA methylation.

DNA methylation in AgRP neurons regulates voluntary exercise behavior in mice.

DNA methylation regulates cell type-specific gene expression. Here, in a transgenic mouse model, we show that deletion of the gene encoding DNA methyltransferase Dnmt3a in hypothalamic AgRP neurons causes a sedentary phenotype characterized by reduced voluntary exercise and increased adiposity. Whole-genome bisulfite sequencing (WGBS) and transcriptional profiling in neuronal nuclei from the arcuate nucleus of the hypothalamus (ARH) reveal differentially methylated genomic regions and reduced expression of AgRP neuron-associated genes in knockout mice. We use read-level analysis of WGBS data to infer putative ARH neural cell types affected by the knockout, and to localize promoter hypomethylation and increased expression of the growth factor Bmp7 to AgRP neurons, suggesting a role for aberrant TGF-ß signaling in the development of this phenotype. Together, these data demonstrate that DNA methylation in AgRP neurons is required for their normal epigenetic development and neuron-specific gene expression profiles, and regulates voluntary exercise behavior.

Early postnatal overnutrition accelerates aging-associated epigenetic drift in pancreatic islets.

Pancreatic islets of type 2 diabetes patients have altered DNA methylation, contributing to islet dysfunction and the onset of type 2 diabetes. The cause of these epigenetic alterations is largely unknown. We set out to test whether (i) islet DNA methylation would change with aging and (ii) early postnatal overnutrition would persistently alter DNA methylation. We performed genome-scale DNA methylation profiling in islets from postnatally over-nourished (suckled in a small litter) and control male mice at both postnatal day 21 and postnatal day 180. DNA methylation differences were validated using quantitative bisulfite pyrosequencing, and associations with expression were assessed by RT-PCR. We discovered that genomic regions that are hypermethylated in exocrine relative to endocrine pancreas tend to gain methylation in islets during aging (R 2 = 0.33, P < 0.0001). These methylation differences were inversely correlated with mRNA expression of genes relevant to ß cell function [including Rab3b (Ras-related protein Rab-3B), Cacnb3 (voltage-dependent L-type calcium channel subunit 3), Atp2a3 (sarcoplasmic/endoplasmic reticulum calcium ATPase 3) and Ins2 (insulin 2)]. Relative to control, small litter islets showed DNA methylation differences directly after weaning and in adulthood, but few of these were present at both ages. Surprisingly, we found substantial overlap of methylated loci caused by aging and small litter feeding, suggesting that the age-associated gain of DNA methylation happened much earlier in small litter islets than control islets. Our results provide the novel insights that aging-associated DNA methylation increases reflect an epigenetic drift toward the exocrine pancreas epigenome, and that early postnatal overnutrition may accelerate this process.

A genomic atlas of systemic interindividual epigenetic variation in humans.

BACKGROUND: DNA methylation is thought to be an important determinant of human phenotypic variation, but its inherent cell type specificity has impeded progress on this question. At exceptional genomic regions, interindividual variation in DNA methylation occurs systemically. Like genetic variants, systemic interindividual epigenetic variants are stable, can influence phenotype, and can be assessed in any easily biopsiable DNA sample. We describe an unbiased screen for human genomic regions at which interindividual variation in DNA methylation is not tissue-specific. RESULTS: For each of 10 donors from the NIH Genotype-Tissue Expression (GTEx) program, CpG methylation is measured by deep whole-genome bisulfite sequencing of genomic DNA from tissues representing the three germ layer lineages: thyroid (endoderm), heart (mesoderm), and brain (ectoderm). We develop a computational algorithm to identify genomic regions at which interindividual variation in DNA methylation is consistent across all three lineages. This approach identifies 9926 correlated regions of systemic interindividual variation (CoRSIVs). These regions, comprising just 0.1% of the human genome, are inter-correlated over long genomic distances, associated with transposable elements and subtelomeric regions, conserved across diverse human ethnic groups, sensitive to periconceptional environment, and associated with genes implicated in a broad range of human disorders and phenotypes. CoRSIV methylation in one tissue can predict expression of associated genes in other tissues. CONCLUSIONS: In addition to charting a previously unexplored molecular level of human individuality, this atlas of human CoRSIVs provides a resource for future population-based investigations into how interindividual epigenetic variation modulates risk of disease.

Epigenetic supersimilarity of monozygotic twin pairs.

BACKGROUND: Monozygotic twins have long been studied to estimate heritability and explore epigenetic influences on phenotypic variation. The phenotypic and epigenetic similarities of monozygotic twins have been assumed to be largely due to their genetic identity. RESULTS: Here, by analyzing data from a genome-scale study of DNA methylation in monozygotic and dizygotic twins, we identified genomic regions at which the epigenetic similarity of monozygotic twins is substantially greater than can be explained by their genetic identity. This "epigenetic supersimilarity" apparently results from locus-specific establishment of epigenotype prior to embryo cleavage during twinning. Epigenetically supersimilar loci exhibit systemic interindividual epigenetic variation and plasticity to periconceptional environment and are enriched in sub-telomeric regions. In case-control studies nested in a prospective cohort, blood DNA methylation at these loci years before diagnosis is associated with risk of developing several types of cancer. CONCLUSIONS: These results establish a link between early embryonic epigenetic development and adult disease. More broadly, epigenetic supersimilarity is a previously unrecognized phenomenon that may contribute to the phenotypic similarity of monozygotic twins.

Maternal exercise during pregnancy promotes physical activity in adult offspring.

Previous rodent studies have shown that maternal voluntary exercise during pregnancy leads to metabolic changes in adult offspring. We set out to test whether maternal voluntary exercise during pregnancy also induces persistent changes in voluntary physical activity in the offspring. Adult C57BL/6J female mice were randomly assigned to be caged with an unlocked (U) or locked (L) running wheel before and during pregnancy. Maternal running behavior was monitored during pregnancy, and body weight, body composition, food intake, energy expenditure, total cage activity, and running wheel activity were measured in the offspring at various ages. U offspring were slightly heavier at birth, but no group differences in body weight or composition were observed at later ages (when mice were caged without access to running wheels). Consistent with our hypothesis, U offspring were more physically active as adults. This effect was observed earlier in female offspring (at sexual maturation). Remarkably, at 300 d of age, U females achieved greater fat loss in response to a 3-wk voluntary exercise program. Our findings show for the first time that maternal physical activity during pregnancy affects the offspring's lifelong propensity for physical activity and may have important implications for combating the worldwide epidemic of physical inactivity and obesity.-Eclarinal, J. D., Zhu, S., Baker, M. S., Piyarathna, D. B., Coarfa, C., Fiorotto, M. L., Waterland, R. A. Maternal exercise during pregnancy promotes physical activity in adult offspring.

Developmental programming of energy balance regulation: is physical activity more 'programmable' than food intake?

Extensive human and animal model data show that environmental influences during critical periods of prenatal and early postnatal development can cause persistent alterations in energy balance regulation. Although a potentially important factor in the worldwide obesity epidemic, the fundamental mechanisms underlying such developmental programming of energy balance are poorly understood, limiting our ability to intervene. Most studies of developmental programming of energy balance have focused on persistent alterations in the regulation of energy intake; energy expenditure has been relatively underemphasised. In particular, very few studies have evaluated developmental programming of physical activity. The aim of this review is to summarise recent evidence that early environment may have a profound impact on establishment of individual propensity for physical activity. Recently, we characterised two different mouse models of developmental programming of obesity; one models fetal growth restriction followed by catch-up growth, and the other models early postnatal overnutrition. In both studies, we observed alterations in body-weight regulation that persisted to adulthood, but no group differences in food intake. Rather, in both cases, programming of energy balance appeared to be due to persistent alterations in energy expenditure and spontaneous physical activity (SPA). These effects were stronger in female offspring. We are currently exploring the hypothesis that developmental programming of SPA occurs via induced sex-specific alterations in epigenetic regulation in the hypothalamus and other regions of the central nervous system. We will summarise the current progress towards testing this hypothesis. Early environmental influences on establishment of physical activity are likely an important factor in developmental programming of energy balance. Understanding the fundamental underlying mechanisms in appropriate animal models will help determine whether early life interventions may be a practical approach to promote physical activity in man.

CpG methylation differences between neurons and glia are highly conserved from mouse to human.

Understanding epigenetic differences that distinguish neurons and glia is of fundamental importance to the nascent field of neuroepigenetics. A recent study used genome-wide bisulfite sequencing to survey differences in DNA methylation between these two cell types, in both humans and mice. That study minimized the importance of cell type-specific differences in CpG methylation, claiming these are restricted to localized genomic regions, and instead emphasized that widespread and highly conserved differences in non-CpG methylation distinguish neurons and glia. We reanalyzed the data from that study and came to markedly different conclusions. In particular, we found widespread cell type-specific differences in CpG methylation, with a genome-wide tendency for neuronal CpG-hypermethylation punctuated by regions of glia-specific hypermethylation. Alarmingly, our analysis indicated that the majority of genes identified by the primary study as exhibiting cell type-specific CpG methylation differences were misclassified. To verify the accuracy of our analysis, we isolated neuronal and glial DNA from mouse cortex and performed quantitative bisulfite pyrosequencing at nine loci. The pyrosequencing results corroborated our analysis, without exception. Most interestingly, we found that gene-associated neuron vs. glia CpG methylation differences are highly conserved across human and mouse, and are very likely to be functional. In addition to underscoring the importance of independent verification to confirm the conclusions of genome-wide epigenetic analyses, our data indicate that CpG methylation plays a major role in neuroepigenetics, and that the mouse is likely an excellent model in which to study the role of DNA methylation in human neurodevelopment and disease.

Independent genomewide screens identify the tumor suppressor VTRNA2-1 as a human epiallele responsive to periconceptional environment.

BACKGROUND: Interindividual epigenetic variation that occurs systemically must be established prior to gastrulation in the very early embryo and, because it is systemic, can be assessed in easily biopsiable tissues. We employ two independent genome-wide approaches to search for such variants. RESULTS: First, we screen for metastable epialleles by performing genomewide bisulfite sequencing in peripheral blood lymphocyte (PBL) and hair follicle DNA from two Caucasian adults. Second, we conduct a genomewide screen for genomic regions at which PBL DNA methylation is affected by season of conception in rural Gambia. Remarkably, both approaches identify the genomically imprinted VTRNA2-1 as a top environmentally responsive epiallele. We demonstrate systemic and stochastic interindividual variation in DNA methylation at the VTRNA2-1 differentially methylated region in healthy Caucasian and Asian adults and show, in rural Gambians, that periconceptional environment affects offspring VTRNA2-1 epigenotype, which is stable over at least 10 years. This unbiased screen also identifies over 100 additional candidate metastable epialleles, and shows that these are associated with cis genomic features including transposable elements. CONCLUSIONS: The non-coding VTRNA2-1 transcript (also called nc886) is a putative tumor suppressor and modulator of innate immunity. Thus, these data indicating environmentally induced loss of imprinting at VTRNA2-1 constitute a plausible causal pathway linking early embryonic environment, epigenetic alteration, and human disease. More broadly, the list of candidate metastable epialleles provides a resource for future studies of epigenetic variation and human disease.

Maternal nutrition at conception modulates DNA methylation of human metastable epialleles.

In experimental animals, maternal diet during the periconceptional period influences the establishment of DNA methylation at metastable epialleles in the offspring, with permanent phenotypic consequences. Pronounced naturally occurring seasonal differences in the diet of rural Gambian women allowed us to test this in humans. We show that significant seasonal variations in methyl-donor nutrient intake of mothers around the time of conception influence 13 relevant plasma biomarkers. The level of several of these maternal biomarkers predicts increased/decreased methylation at metastable epialleles in DNA extracted from lymphocytes and hair follicles in infants postnatally. Our results demonstrate that maternal nutritional status during early pregnancy causes persistent and systemic epigenetic changes at human metastable epialleles.

Major epigenetic development distinguishing neuronal and non-neuronal cells occurs postnatally in the murine hypothalamus.

Prenatal and early postnatal environment can persistently alter one's risk of obesity. Environmental effects on hypothalamic developmental epigenetics constitute a likely mechanism underlying such 'developmental programming' of energy balance regulation. To advance our understanding of these processes, it is essential to develop approaches to disentangle the cellular and regional heterogeneity of hypothalamic developmental epigenetics. We therefore performed genome-scale DNA methylation profiling in hypothalamic neurons and non-neuronal cells at postnatal day 0 (P0) and P21 and found, surprisingly, that most of the DNA methylation differences distinguishing these two cell types are established postnatally. In particular, neuron-specific increases in DNA methylation occurred extensively at genes involved in neuronal development. Quantitative bisulfite pyrosequencing verified our methylation profiling results in all 15 regions examined, and expression differences were associated with DNA methylation at several genes. We also identified extensive methylation differences between the arcuate (ARH) and paraventricular nucleus of the hypothalamus (PVH). Integrating these two data sets showed that genomic regions with PVH versus ARH differential methylation strongly overlap with those undergoing neuron-specific increases from P0 to P21, suggesting that these developmental changes occur preferentially in either the ARH or PVH. In particular, neuron-specific methylation increases at the 3' end of Shh localized to the ARH and were positively associated with gene expression. Our data indicate a key role for DNA methylation in establishing the gene expression potential of diverse hypothalamic cell types, and provide the novel insight that early postnatal life is a critical period for cell type-specific epigenetic development in the murine hypothalamus.

Early postnatal nutrition determines adult physical activity and energy expenditure in female mice.

Decades of research in rodent models has shown that early postnatal overnutrition induces excess adiposity and other components of metabolic syndrome that persist into adulthood. The specific biologic mechanisms explaining the persistence of these effects, however, remain unknown. On postnatal day 1 (P1), mice were fostered in control (C) or small litters (SL). SL mice had increased body weight and adiposity at weaning (P21), which persisted to adulthood (P180). Detailed metabolic studies indicated that female adult SL mice have decreased physical activity and energy expenditure but not increased food intake. Genome-scale DNA methylation profiling identified extensive changes in hypothalamic DNA methylation during the suckling period, suggesting that it is a critical period for developmental epigenetics in the mouse hypothalamus. Indeed, SL mice exhibited subtle and sex-specific changes in hypothalamic DNA methylation that persisted from early life to adulthood, providing a potential mechanistic basis for the sustained physiological effects. Expression profiling in adult hypothalamus likewise provided evidence of widespread sex-specific alterations in gene expression. Together, our data indicate that early postnatal overnutrition leads to a reduction in spontaneous physical activity and energy expenditure in females and suggest that early postnatal life is a critical period during which nutrition can affect hypothalamic developmental epigenetics.